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(~| . Abstract 

"y ' The aim of this paper is to derive a raw Bloch model for the interac- 

C^ . tion of light with quantum boxes in the framework of a two-electron- 

I i ' species (conduction and valence) description. This requires a good 

understanding of the one-species case and of the treatment of level 
^^ ' degeneracy. In contrast with some existing literature we obtain a Li- 

. J^ . ouville equation which induces the positiveness and the boundcdncss of 

QQ ' solutions, that are necessary for future mathematical studies involving 

(yQ . higher order phenomena. 

p; 

^ ■ 1 Introduction 

oo : 

^D ■ Bloch equations are a very common model to describe the time evolution of 

^ , a system of electrons in different contexts such as gas of electrons, glasses, 

'k^ I crystals. Mathematical results have been obtained in these contexts which 

^ ■ ensure in particular that the variable of Bloch equations, namely the density 



matrix, keeps through the time evolution some properties, such as hermicity 
and positiveness (of diagonal elements and of the matrix as an operator) [H 
[2]. In the above contexts, the indices in the density matrix are integers which 
distinguish between the different electron levels. Diagonal entries model the 
population of the levels, while off-diagonal entries model coherences between 
two levels. 

Bloch equations have also been derived for quantum wells [5], [6] . They 
look very similar but the variables are different. In the quantum well con- 
text, we distinguish between two species of electrons (valence and conduc- 
tion electrons) and variables are indexed by (two-dimensional) wave vectors. 



Electrons can interact directly only if they have the same wave vector, which 
induces intra-band coherences to be zero and valence and conduction elec- 
trons to be coupled two-by-two. This is not a mathematical problem since 
at first order, i.e. without considering relaxation effects, such a system is a 
juxtaposition of two-level models. 

The model for quantum wells has been extended to quantum boxes [1]. 
Due to the three-dimensional confinement of electrons, variables are now 
again indexed by integers. Valence and conduction electron levels are not 
coupled two-by-two any more and any valence electron can interact with any 
conduction electron. But in the propounded model intra-band coherences 
are still considered as zero. This destroys the Liouville structure of the 
system which in other respects resembles much the gas model. Another 
puzzling point is the fact that levels are supposed to have any degree of 
degeneracy. This is not much coherent with the Pauli exception principle 
that electrons (which are fermions) are supposed to follow. This induces 
also mathematical and modelling problems which we also address here. 

In this paper we are interested in the derivation of a raw Bloch model 
for a quantum box system. This model will be derived carefully from com- 
putations involving creation and annihilation operators for both species of 
electrons. Our goal is to obtain a model with appropriate mathematical 
properties in order to be able to prove the conservation of certain physical 
properties through time. Our study will be restricted to the leading order 
terms in the Hamiltonian, leaving extra contributions such as Coulomb ef- 
fect or electron-phonon interaction to future works. More precisely this is 
not a restricted study but the first step towards a broader study, since the 
properties proved on the first order terms are necessary to study the higher 
terms. 

The outline of this paper is as follows. In Section [2] we derive carefully 
a one-species Bloch model such as the one for gazes. The structure of this 
model is analysed. Positiveness results are recalled and boundedness of the 
populations is discussed. This motivates Section [3] where the problem of 
degeneracy is thoroughly addressed. Section U] is devoted to our central 
problem, i.e. the derivation of a two-species Bloch model. This model has 
the desired mathematical properties. The details of the computations and 
an analogy with a boson model are postponed to various appendices. 



2 The case of one species of electrons 

2.1 Commutators and Heisenberg equation 

Let A and B be two operators, we define their commutator by 

[A, B]= AB- BA 

and their skew-commutator by 

{A,B} = AB + BA. 

For an operator A, we will define the associated observable {A) 

{A) = Tv{SoA) 

by averaging with respect to the initial state density Sq of the system. If 
the system is described by an Hamiltonian H, the time-evolution for this 
observable is given by the Heisenberg equation 

ihdt{A) = {[A,H]). 

Bloch equations are the Heisenberg equation of motion when the observ- 
able is the density matrix. 

2.2 Operators and commutation rules 

We suppose there is only one species of electrons, like in the case of gas 
atoms, and that levels are indexed by integers i G /. For the ith level, 
we define the creation and annihilation operators cj and a. Electrons are 
fermions and should respect the Pauli exception principle. The correspond- 
ing skew-commutation rules are 

{ci, Cj} = {cl,clj} = 0, {q, c]} = 6ij. 

This implies in particular that 

this last equality meaning clearly that it is impossible to create twice or 
annihilate twice the same electron. This is the Pauli exception rule. Let A 
and B be the products of ua and ns such operators which all skew-commute 
two-by-two, then 

J 2AB if UA or ns is odd, 
^^^""^ = [0 else. 

Since we only use two-operator observables, we will mostly always be in the 
second case. 



2.3 Hamiltonian and observables 

We define the Hamiltonian as the sum of a part due to the electrons and 
only, and a part due to their interaction with an external electric field E(i): 
H = H'' + H^ with 

kei 
where e^ is the energy of an electron in the kth level and 

{k,l)€P 

where M^; is an element of the dipolar moment matrix M. The 1/2 coeffi- 
cient is due to the fact that interactions are counted twice in this sum. The 
dipolar moment M^^ can be expressed as (V';*|r|V'A:) in terms of the wave 
functions associated to each level and the local position vector r. This im- 
plies in particular that the case k = I does not contribute to the sum since 
Myfcfc = 0. We will also use the fact that M is hermitic : M^^ = M^^. 

The observable we are interested in is the density matrix, which ele- 
ments are the pij = {CjCi). It is clear that this matrix is Hermitian. In the 
computations, we distinguish between diagonal terms of this matrix also 
called populations pu = (clci) and off-diagonal terms, pij, i j^ j also called 
coherences. 

2.4 Computation of the raw Bloch equations 

We have to compute the commutators of the density matrix with both 
Hamiltonians. The details of the computations are given in Appendix |Al 
We obtain 

{[c]ci,H'']) = {ei -ej)pij, 

and 

2{[c]ci,H^]) = ReE(t) • Y^i^ikPkj - MkjPik)- 

kei 

Gathering these results and the Heisenberg equation, raw Bloch equations 
read 

ihdtpij = {ei - ej)pij + ReE(t) • ^(MjfcPfcj - M^k^Pik)- 

kei 

If we denote E = diag({ej}ig/), this can be cast as a Liouville equation 

[fidtp=[V{t),p], wherey(t) =^ + ReE(t)-M. (1) 



In Appendix[B]we show that considering electrons as bosons leads to exactly 
the same equation, although the intermediate computations are different. 
This is not true any more if you consider extra contributions to the Hamil- 
tonian such as Coulomb effect or electron-phonon interaction, but this is an 
other story. 

2.5 Positiveness and trace results 

Equation ([T]) clearly preserves the Hermitian structure of p. We can give an 
exact solution to equation ([1]), namely 

p{t) = exp {-'- ^ V{T)dT^ p(0) exp {^^ ^ V{T)dT^ . 

This evolution equation preserves positiveness. Indeed, if /j(0) is a non 
negative operator, that is for all X G /^(/), X* p{Q)X > 0, then 

X*p{t)X = X*exp (-^ I V{T)dT\ p(0)exp (^ f F(r)dT j X 

= Y*p{0)Y, where y = exp (^ / F(r)dT J X 

Therefore X*p{t)X > 0. A corollary of this is that for all time t, and all 
level j, pjj{t) > 0. 

The total population is the sum of populations, which can be described 
by the trace of the density matrix. Since commutation is a trace preserving 
operation, we have 

indtTTp = Tt[V{t),p] =0. 

Hence Trp(t) = Ti p{0) = A'tot- The raw Bloch equations are trace pre- 
serving. In other words the total number of electrons is preserved through 
the time evolution. This together with non negativeness ensures that for all 
time t 

< pjj{t) < Ntof 

We want also to know whether each level remains bounded by 1, which 
would be the expression of the Pauli exception principle. This is indeed the 
case. If we set p = Id — p, where Id is the identity matrix with the same 
dimensions as p, then 

[VM-p] = [VM]-[v,p] = -[v,pi 

dt(ld-p) = -dtp. 

Therefore p is also solution to the raw Bloch equations, and in particular 
Piiit) ^ 0, which means pii{t) < 1. 



3 Degenerate levels 

In [3], levels are supposed to be degenerate, with a degeneracy order that 
depends on the levels. In a fermion description levels cannot be degenerate 
since each level population is bounded by 1. This would a priori be possible 
in a boson description. Let us first investigate this situation. 

The same proof as above shows that if p is solution to the raw Bloch 
equations and the level populations are initially bounded by the same con- 
stant (say d), then pii{t) < d. To derive a d-degenerate raw Bloch model, it 
is therefore possible to consider electrons as bosons and impose a condition 
on the initial data. But this solution is not very satisfactory from two points 
of view. First this cannot be extended to the case when the degeneracy de- 
pends on the level, i.e. the ith level is dj-degenerate. Second if we want 
to look at other contributions in the Hamiltonian, bosonic and fermionic 
commutation rules do not lead to the same model. 

The question is therefore: can we derive a degenerate model from fermion 
calculations with level dependent degeneracies? 

3.1 Deriving a degenerate model 

For z G / and n = l,...,di, we denote by c" and cf the creation and 
annihilation operators associated to the nth degenerate sub-level of level 
i. To ensure a coherent degenerate model we have to assume the following 
commutation rules 



"t ^Mi _ n /^" ^'"ti 



{c7, cf} = {c-\cf'} = 0, {cr, cf'} = 6,,5n 



rri' 



The energy ej does not depend on the sub-level index n. Since the levels 
are exactly the same (and not just have the same energy), the associated 
wave functions are the same and therefore the dipolar moment entries only 
depend on the levels and not on the sub-level. We still denote them by Mjj. 
The Bloch equations governing pff' = {cf cf) is the Full Degenerate Bloch 
(FDB) model 

d-k 
fcG/ P=l 

We notice that natural quantities that arise in the right handside are 

Pkj — z_^ Pkj ' Pik — z_^ Pik ' 

p=i p=i 



and we write the equation for p" 



«j 



mtpf = ihdt^p: 



nm 
m=l 



m=l k£l m=l m=l 

= (e, - e,)p^; + ReE(t) • E(M.fePt+ - Mk,d,pl+), 

k£l 

where we have introduced 

di dj 
n=l m=l 

and we obtain a closed set of equations for the density matrix p'^^, which 
is the Condensed Degenerate Bloch (CDB) model 

n=l 

= (e. - e.) E pl^ + i^«E(t) • E(M.. E pt^ - M/c.rf. E pi''^ 

n=l kal n=l n=l 

= (e, - e,)p++ + ReE(t) • J](M,fcd,p++ - Mkjdjp++). 

k£l 

It is not exactly the same equation as for non degenerate levels since the 
degeneracies occur in the equation coefficients. If we set 

pr 



aij = -j===, Njj = Mijy/didj, 

we recover the usual Bloch equation 

ihdtaij = {ei - ej)aij + Re E(t) • E(NifcO-fcj - ^kj(^ik)- (3) 

kei 

This ensures in particular that a (and therefore p^~^) defines a non negative 
operator if this is valid at the initial time. 



3.2 Boundedness of degenerate levels 

The diagonal elements are an = pfi^ /di and we would expect this model to 
ensure aii{t) < 1 (i.e. p^^ < di) for all time. By the same arguments as 
above applied on equation ([3]), this is true if this hold at the initial time. 
The problem is that this condition is not natural when dealing with the 
variables of the FDB model ^. Indeed the diagonal entries of the CDB 
model are not the sum of the populations of a given level but also include 
intra- level coherences. An other consequence of the non negativeness of p is 
that \pij\ < ^PiiPjj (for the variables of equation ([1]), see [2]). If plugged in 
the definition of an this yields 




di , , 

n=l m=l 



rii 



m=l / 



Since there can be configurations where an = di, the only way to ensure that 
(^iiiO) < 1 is to impose that intra-level coherences are initially zero. This 
certainly is valid in most experimental situations. Under this vanishing 
condition, the CDB model preserves the natural property p^~^ < di. In the 
opposite situation, it is not possible to preserve the property p^^ < d? which 
reverts to preserve an < di and is the same problem as deriving a degenerate 
model based on a boson derivation of Bloch equations, the initial problem 
that we have eluded. 

4 The case of two species of electrons 

4.1 Conduction— valence notations 

To be closer to the above computations, we will first present the model of 
two species of fermions considering electron operators (and not electron and 
hole operators, as in ^). For i € P (a set of integers, indexing conduction 
electrons), we keep the notations cj and Cj for the creation and annihilation 
operators associated to conduction electrons. For i £ I^ (indexing valence 
electrons), we will use the non conventional notation vj and Vi for the cre- 
ation and annihilation operators associated to valence electrons. Along with 
the already defined commutation rules for conduction electrons 



{cj,Cj} = {c|,c|} = 0, {ci,c\} = 6., 



«j' 



we have the same rules for valence electrons 

{vi,Vj} = {vl,vj} = 0, {vi,vj} = 6ij 

and commutation rules between the two species 

[ci,Vj] = [cl,vj] = [Ci,vj] = [cl,Vj] = 0. 
We now consider the Hamiltonians 

H"^" = I E iW)-Mlivlvi+E*{t)-Ml^v}vk), 
^Lcv ^ ^ (E(t).M|l4t;; + E*(t).M^r^/cfc)- 



The density matrices we are interested in are intra-band densities 
and inter-band densities 



p'ij = {v]c,), Pi;i = {c]vi). 

We will not write the equations for p™ but an appropriate definition of M^ 
ensures that 

P P 

is Hermitian. 

We now have to compute the commutators of these matrices with the 
three above Hamiltonians. This is done in Appendix [Cl 



4.2 Towards a two-species Bloch equation 

Gathering all the computations of Appendix O we obtain 

+ E(i) • Y. mTkPl) - E*(t) • Y ^TkPTk. 

kel'' kGl^ 



mdtpl = (el - e])pl + Re E(t) • ^ (MJ.p^^. - MI^-pJ,) 

-E(t).5^M^,>r,= -E*(t).5]M: 

ihdtp^ = {el - e])pT, + Re E(t) • ( ^ M^,p^,] - ^ ^l^PTk 



cv* cv 
ki Pkji 



cv 

V / ) 



If we denote E'^ = diag({ej}jg/c), E"" = diagdej'jjg/v) and 

/ -E^ + ReE(t)-M^ E(t)-M' 

»/W-\^ E*(t)-M'=" S^ + ReE(t) -M 

we have the Bloch equation in Liouville form 

indtp'°' = [v{t),p'°\ 

We can in a straightforward way apply our discussion about the degeneracy 
of levels to this two-species context. 

4.3 Electron— hole formulation 

Since all valence electrons do not play a role in the interaction with con- 
duction electrons, only those who interact are described and not by their 
presence but their absence: a hole in the valence band. We formulate anew 
the former notations and results, denoting for i G / = /^ by dj = Vi the 
creation operator of a hole (annihilation of the corresponding electron) and 
di = v\ its annihilation operator. The hole energy is the opposite of energy 
of the corresponding electron e^ = —ej. Moreover d\di = l — did] = l — vjvi, 
hence the Hamiltonian for holes is derived from the Hamiltonian for valence 
electrons via 

^^ = E '^I'^k = E -^'(1 - 44) = E -^' + E ^'44. 

k&I" k£l^ ki^I^ k£l^ 

10 



The first term is a constant that will not play any role in the commutators, 
hence we define 

The intra-band and inter-band density matrices involving holes are defined 
by p}j = {djdi) = 5ij - pjj and p^j = (djCi) = p'fj . The system for these 
variables read 



4) 



ihdtplj = (e? - ^'j)p'io + ReE(t) • Y, i^lkPlj - ^IjPl 

kei'= 

+ E(t) . ^ Mif^ - E*(t) . J2 MffcV^L 
ihdtp^^ = {e^ - e^)p'^^ + ReE(t) • J^ i^]kPl " M^p^,) 



(4) 



„ch\ 



ihdtp^ = (ef + e^)f^ + ReE(t) • (^ M^,/,g - ^ Ml^.p?,^; 

+ E(t) • ( ^ M^,^<5,, - p%) - Y ^jf^k)- 

k£l^ kel" 

which has lost the Liouville structure. Here we have chosen to denote 
MJI- = MJ- and Mjj^ = M?J but no other choice would help to recover 
this structure. Mathematical results have to be obtained from the electron- 
electron formulation. The electron-hole formulation should be only kept 
for simulation or intuition in introducing new terms such as electron-hole 
recombination. 

4.4 Vanishing intra-band coherences 

In comparison with the model proposed in [4], a natural question is now 
whether we recover their model by imposing vanishing intra-band coher- 
ences, i.e. /7?- = pjj = for z 7^ J, . In this reference the variables are 

"■I ~ Pii^ ^^ ~ P^j ^^^ Pji ~ Pfj (with a subtlety about coupling with 
forward and backward propagating optical fields, which we do not separate 
here). With this set of variables and vanishing intra-band coherences, the 
system ([4]) now reads 

mtnl = E(t). J^Mf,hp^,-E*(t). J]Mf,>fc„ 

11 



ihdtpji = (ef + e^)pji + Re E(t) • ( J^ M^^p^-fc - ^ M^^pfc,) 

The difference witli tlie equations announced in [1] (apart from extra terms 
describing phenomena out of the scope of the present paper) are the terms 
involving M'^ and M . 

Instead of assuming that intra-band coherences are zero, we could as- 
sume that M*^ and M are zero (or, since this is certainly not true, that the 
Hamiltonians H^'^ and H^^ should not be taken into account in the deriva- 
tion), we also then do not recover the equations in ^ because of extra terms 
in the equation for pji. Under this only assumption intra-band coherences 
are not zero, even if taken to be zero at the initial time. 

The equations in ^ have therefore lost the Liouville structure and are 
not suitable for the mathematical analysis. The model we propose is far 
better in this respect. If assuming zero intra-band coherences is not coherent 
with the Liouville structure, we can easily have them very small by simply 
assuming them to be initially zero and adding in their equations strong 
damping terms. As shown in ^, positiveness results can still be obtained 
in this context. Besides the frequencies of E are chosen to be close to the 
gap frequency, in order to match valence-conduction transition frequencies. 
The intra-band coherences are hence not directly excited by such a wave. 

5 Conclusion 

We have derived raw Bloch equations for quantum boxes, i.e. two species 
of electrons. These equations have the Liouville form which allows to prove 
positiveness results. We have also discussed the problem of level degen- 
eracy and derived a condensed degenerate Bloch model in which diagonal 
entries remain bounded by the corresponding level degeneracy through time 
evolution, even when the degeneracy order is level-dependent. 

Positiveness and boundedness results are necessary to prove similar prop- 
erties when extra terms are added in these raw equations. They are also 
required to prove sharp existence results when coupled with a model for 
wave propagation (e.g. Maxwell equations) as already done for one species 
of electrons ^ . 

Future work will consist in introducing with the same care extra con- 



12 



tributions to the Hamiltonian to model Coulomb effect, electron-phonon 
interactions, electron-hole recombination, . . . Our goal is to obtain a model 
close to that of [3] but with anew possible other contributions, and full 
equations for intra-band coherences. 

A Detail of one species fermion computations 

Let us first compute the commutator with the free electron Hamiltonian. 
We compute separately the cases i i^ j and i = j but notice afterwards that 
results reads the same in both cases namely 

([cJci,iJ'^]) = ^efc(cjci4cfc-4cfcc{ci) 

k 

= ei{c\cic\ci - c\cic\ci) = 0, 

= ei(c]cj) -ej(c]Q) 
= (ei - ^j)pij- 

Now we compute the commutators with H . 

2{[c\c,,H^]) = E(t) • Y, ^iMc4ci - cJq4q) 

+E*(t) • ^M*;(4Q4ci - c\cic\ci) 

+E(t) • ^ Mfci(c|cicI,Ci - (^^Ci<^iCi) 

ki^i 

+E*(t) -Yj^UcWVk - c\ckc\ci) 

k^i 

= E(t) ■Y.'^iMci) - E*(t) • J;M*,(4c,) 

-E{t) . Y, M,,{clc,) + E*(t) • Y m^{4ck) 

k^i ky^i 

([cJci/L]) = ReB{t)■Y^^kPk^-ReE{t)■Y^ikPik 

k k 
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2([c]q, i^L]) = E(t) • Y^ Ma{c]c4ci - c\cic]c^) 

+E*(t) ■Y,^]i{c]cic\cj - c\cjc]ci) 
+E(t) • ^ Mfcj(c]ci4cj - 4^4^) 

+E*(t) • j]M^,(4Q4cfc - 4cfc4Q) 
+E(t) • Mjj(4ci4cj - 4^4^^*^ 

+E*(t) • M*,(4q4c, - c\cjc\c,) 

= E(t) • ^ Mii(4(l - 4Q)q - 4^/4^*) 

+E*(t) • J;M*,(4q4c, - 4(1 - c\cj)c,) 

+E(t) • ^ Mfcj(4Q4ci - 4(1 - 4'^i)''^) 

+E*(t) • ^M^,(4(l - 4Q)cfc - c\cJ^o>^ 

+E(t) • Mij-(4(1 - c\ci)cj - 4(1 - c]cj)c^) 
+E*(t) • M*,(4(l - 4ci)c, - 4(1 - 4c,)q) 
= E(t) • J] M,i(4Q) - E*(t) • J] M*,(4q) 

-E(t) • ^ Mfe,(4c,) + E*(t) • j; MUclck) 

k+i kytj 

+E(t) • MijicU - 4q) + E*(t) • MUcUj - 4c, 



([4q, H^]) = ReE(t) • ^ M,kPkj - ReE(t) • ^ M^jP, 






fc 



B Boson computations 

If electrons are considered as bosons, or in other words without Pauh excep- 
tion rule, the definitions for the operators, Hamiltonian and density matrix 
are the same as for fermions, but the computations differ because of com- 
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mutation rules between creation and annihilation operators which are now 

[ci,Cj] = [c], c]] = 0, [a, c]] = 5ij. 

We still have {[clci, H'^]) = 0. We first compute the commutators of the 
density matrix elements with the electron Hamiltonian 

{[c]ci,H'']) = y^ ek{c\cic\ck - c\ckc]ci) 

k 
~~ ^ix^j^i^i^ C^^CiCjCij + Cj\CjCiCjCj CjCjCjCi) 

= ej(c](l + c|ci)Q - c|Qc]ci) + ej{c\cic\cj - c](l + cyj)ci) 
= eiic^i) -ejic^jCi) 

Then we compute the commutators of the density matrix elements with the 
laser Hamiltonian 

2([c1q, i/L]) = E(t) • Y, ^^Mc4cl - c\cic\c,) 



+E*(t) • J]M*,(4q4q - c\c4c. 



l^i 



+E(t) • ^ Mfcj(c|cic|.Cj - 4cic|ci 



k^i 



+E*(t) .Y,m^{cU^cUk - c\ckc\c^) 



k^i 



E(t) • ^M,i(cJ(l + c\ci)ci - c\cic\ci^ 



i^i 



+E*(t) ■Y.micUicU - cj(l + 4cOc* 



i^i 



+E(t) • Y^ Mkiiclacla - 4(1 + c|q)q 



k^i 



+E*(t) • ^M^,(cJ(l + 4c,)ck - c\ckc\ci^ 



k^i 



* iJ^ 



E(t).^M,,(4q)-E*(t).^M*,(c 
-E(t) .Y,Mk^{cU)+B*{t) • J^M^,(c|cfc) 

k^i ky^i 
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E(t) • ^ mikPki - E*(t) • ^ mkiPik 

-E{t) ■ Y, MkiPik + E*(t) • Y, ^ikPki, 

ky^i kj^i 

ReE(t) ■Y^^kPk^ - ReE(t) -^^IkPik, 

k k 

E(t) • Y Mii{c]ciclci - c\cic]ci) 

+E(t) • Y ^kj{c]ciclcj - c\cjc]ci) 

+E(t) • Mij(c]cic|cj - cicjcjci) 
+E*(t) • M*,(4qcJc, - cJc,4q) 
E(t) • ^Mij(c](l + cJcOq - cJqc]q) 

+E*(t) • ^M*,(4q4c, - c|(l + 4c,)q) 

+E(t) • Y ^kj{c]ciclcj - 4(1 + cjcjOci) 

ky^i 

+E*(t) • J;M^,(c](1 + clQ)cfc - 4ckc]c,) 

+E(t) • Mij{c]{l + 4ci)cj - 4(1 + c]cj)ci) 
+E*(t) • M;(4(1 + 4c,)c, - 4(1 + c]c,)c,) 
Bit) • J]M,,(4q) - E*(t) • J]M*,(4q) 

-E(t) • j;Mfc,(4Q) +E*(t) -YMU^^k) 

ky^i kytj 

+E(t) • Mi,-(c]cj - 4q) + E*(t) • M*i{c]cj - c\c^ 
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= E(t) • ^ MikPkj - E* (t) • Y, ^kjpik 

-E(t) • Y, Mfciftfc + E*(t) • Y ^rkPkj 
k^i kytj 

+E(t) • Mij{pjj - pii) + E*(t) • M.ij{pjj - pii), 
{[c]ci,H^]) = Re E(t) • Y ^ikPkj - Re E(t) • ^ MkJP^k■ 

k k 

Although calculations are different, we find the same results as in the 
boson case. If we consider the Pauli exception rule or not we obtain the 
same raw Bloch equations. 

C Detail of two-species fermion computations 

Because conduction and valence operators commute, we clearly have 

{[c]c,,H^]) = {[v}v,,H^]) = {[c]c,,H^^]) = {[v}v,,H^^]) = 0. 
Besides we have already computed in the one-species case 

([c]q, H'^]) = {el - e5)p?,, {[v]v,, W]) = {ej - e])pl, 

{[c]c,,H^'']) = Re^{t)-Y{M.'i,pl^-ml^p1,), 

kel<= 
k&I" 

We now compute 

{[v]ci, H"]) = Y 4(^Jci4cfe - cl(^kv]ci) 
kei'= 

= 4{v]cic\ci - c\civ]ci) = e1{v]ci) = e^ipf], 
{[v^a, H"]) = Y 4{v]civlvk - vlvkvyi) 

k£l^ 

v/tt t t\ v/t\ vcv 

= ej {v'jdv'jVj - v'jVjv'jCi) = -ej (i;jq) = -e^p^j . 

These computations are valid also for i = j, which is a case which has to be 
treated apart for the following computations. 

(k,i)ei''xi^ 

+E*{t) . Yl M.f^*{4c,vjck - vjck4c,) 
{k,i)ei''xr 

17 



= E(t) • J2 ^7{4vi) - E*(t) • J2 ^Tivjc,) 
1(^1^ lei" 

{[c]c,,H'^n) = E(t). J2 Mti{c]c4vi-clvic]c,) 

{fc,/)e/<:x/^ 

+E*(t)- Yl Ml^i^clc^vjck-vjckcla) 
{k,i)ei'=xi'' 

= Bit).Y,^u{4c^4vl-4vlc]c,) 
lei" 



+E(i) • Y, M-^Jiclciclvi - c]vic]ci) 
lei" 

+E*(t) • Y ^Ticlcivjc^ - v\cic\c,) 
ler 



lei" 



E(t) • Y ^li^vi) - E*(t) • Y M,T( 



lei" lei" 

\ . 

3 



V'lCi) 



E(t) • Y ^uPh - E*(i) • Y ^TPu 



lei" lei" 

{[vU^H"^'^]) = E(t). Y ^H{v\v^4vl-4yl4v^) 

{k,i)ei'=xi" 

+E*(t). Y M.'j^mvlv.vjck-vjckvlv, 
{k,i)ei'=xi" 

= E(t) • Y ^tM^^4^^ - 4v.4u:) 

kei" 
+E*(t) • Y Mr(4^'.4cfc - v\ckvU) 



kei'' 
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-E{t) . Y, ^tipjk + E*(t) • Y. M^fciVfci 






+E*(t) • Y ^fi*^^h^W - v\ckv]v^) 
{k,i)ei<=xi'' 

Bit).Y^T^{v}y^cU-clv,v]v,) 

kei'= 



+E(i) • j; M^,](t;]z;,4z;, - clv,v]t 

k&I'^ 

+E*(t) • Y ^t:{vhv\ck - v\cuv]v,) 
kei" 

+E*(t) • Y ^fi'{v]viv]ck - v]ckv]vi) 
kei'= 

-E(i) • Y ^TjicU) + E*(t) • Y ^ti*{v}ck) 
kei" kei<= 

-E(t) • Y ^kjPTk + E*(i) • Y M'^r/^fc}- 



kei'^ kei" 

Once more by different computations we obtain the same result for both 
cases i = j and i j^ j. Now 

2{[v]c,,H'^']) = E(t). Y MUv]c4ci-ciciv}c,) 

(fc,0e(/<=)2 

+E*(t) • ^ Ml*i{v]ci4ck - c\ckv]ci) 
{k,i)e{i<=)2 

= E(t) • Y ^a{v}c4ci - c\civ]o^ 



+E(t) • Y ^lAvhA<^^ - clciv]ci) 
kel'= 

+E*(t) • Y M-r(fjQc|"ci - 4civ]ci 



lei" 



+E*(t) • Y Ml* {v}c^4ck - c\ckv]ci 



k&I" 
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= 2 Re E(t) • ^ M^,(^Jq) = 2 Re E(t) • J] M^,p=,], 



lei 

2{[v}c„H^^]) = E(t). Yl mi{v]c,vlvi-vlviv} 

(fe,0e(/-)2 



+E*(t) • Yl M[*0e(/v)2(^Jc.^|- - -^--^ 



(fc,0e(/^)2 
E(t) • ^ M]i{v]civ]vi - v\viv] 



vlviv'jCi) 

Vk - vjvkvjci) 



lei'' 



J.,. J, 



+E(t) • Y MIj^JQ^t;, - vlvjvja) 

k£l^ 

+E*(i) • Y ^]t{v]c4v, - v\vjv]c,) 
i£i^ 

+E*(t) • Y ^TMciv]vk - v]vkv]ci 



-2ReE(f).^MI^.(4Q) = - 
kei^ 



2ReE(t)- ^M; 

fcG/- 



■V cv 
-kj rik • 



Last 



+E*(t)- Y. ^Ti*{vhv\ck-v\ckv]ci) 



{k,l)el'=xl^ 

E{t).Y^i{^y4^i-4viv] 



Vlv!jCi) 



lei^ 
ker- 



u 



-B{t).ml]{v]c4v,-c\v,v]c,) 

+E*(t) • Y ^T{vhv]ci - v\c^v]ci) 
lei-' 

+E*(t) • Y Mfc7(i']c.«Jcfe - v]ckv]ci) 
kei" 

-E*(t) • M'r]*{v]c^v}c^ - v]c^v]c^) 
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B{t)■Y,M'r[{v]vl-6,lclc^) 

lei" 



kei'= 
-E{t).MrVvlv^-cU 



Bit) ■ Y^ M'ilivjvi) - Bit) . Y^ Mf^icU) 
lei'' kei''- 



kei" kei" 
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